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Abstract. - The behavior of two unidirectionally coupled chaotic oscillators near the gen- 
eralized synchronization onset has been considered. The character of the boundaries of the 
generalized synchronization regime has been explained by means of the modified system ap- 
proach. 



Chaotic synchronization is one of the fundamental nonlinear phenomena actively studied 
recently [1]. Several different types of chaotic synchronization of coupled oscillators, i.e. gen- 
eralized synchronization (GS) [2], phase synchronization (PS) [1], lag synchronization (LS) [3], 
complete synchronization (CS) [4] and time scale synchronization (TSS) [5] are well known. 
There are also attempts to find unifying framework for chaotic synchronization of coupled 
dynamical systems [6-9]. 

One of the interesting and intricate types of the synchronous behavior of unidirectionally 
coupled chaotic oscillators is the generalized synchronization. The presence of GS between the 
response x r (t) and drive x^(f) chaotic systems means that there is some functional relation 
x r (t) = F[xd(f)] between system states after the transient finished. This functional relation 
F[-] may be smooth or fractal. According to the properties of this relation, GS may be divided 
into the strong synchronization and weak synchronization, respectively [10]. There are several 
methods to detect the presence of GS between chaotic oscillators, such as the auxiliary system 
approach [11] or the method of nearest neighbors [2,12]. It is also possible to calculate the 
conditional Lyapunov exponents (CLEs) [10, 13] to detect GS. The regimes of LS and CS are 
also the particular cases of GS. 

One of the interesting aspects of the generalized synchronization study is the analysis of 
the onset of this regime. In particular, the intermittent behavior is known to be revealed on 
the onset of GS [14,15] as well as in case of the lag [3,16,17] or phase [18-20] synchronization. 

At the same time, there are known examples of the unidirectionally coupled chaotic systems 
for which the location of the generalized synchronization onset on the "parameter mismatch 
- coupling strength" plane differs radically from the other synchronization types. Indeed, 
for two unidirectionally coupled Rossler oscillators with identical control parameters the value 
of the coupling strength corresponding to the onset of GS is twice as much as for the same 
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oscillators with parameters detuned sufficiently [21]. For two unidirectionally coupled one- 
dimensional complex Ginzburg-Landau equations with sufficiently detuned control parameters 
the threshold of the GS regime onset does not depend on the value of the drive system 
parameter when the response system parameters are fixed [22]. Alternatively, for all other 
types of chaotic synchronization the dependence of the threshold of the synchronous regime 
arising on the value of the control parameter mismatch behaves in the different way, i.e. when 
the control parameter mismatch decreases the coupling parameter value corresponding to the 
onset of the synchronous regime also declines. 

This letter aims to explain the mechanisms resulting in the generalized synchronization 
arising in the unidirectionally coupled chaotic oscillators and determining the location of the 
GS onset. 

The causes of the generalized synchronization arising may be clarified by means of a mod- 
ified system approach proposed in [23] . Let us consider the behavior of two unidirectionally 
coupled chaotic oscillators with slightly mismatched parameters 

±d = H(x d ,g rf ) 

x r = H(x r ,g r ) +eA(x d - x r ), 

where the state vectors of the drive and response systems, respectively; H defines 

the vector field of the considered systems, gd and g r are parameters vectors, A = {Sij} is 
a coupling matrix (6a = or 1, 8{j — (i ^ j)), e is a scale parameter characterizing the 
coupling strength. 

In this case one can see that the response system x r (t) may be considered as a modified 
system 

x m (i)=H'(x m (t), E ) (2) 
(where H'(x(t)) = H(x(i)) — eAx(i)) under the external force eAxd(t): 

u m (i) = H'(u m (i),e) +eAx d (i), (3) 

It is easy to see that the term — eAx(t) brings the additional dissipation into the system J2J. 
Indeed, the phase flow contraction is characterized by means of the vector field divergence. 
Obviously, the vector field divergences of the modified and the response systems are related 
with each other as 

N 

divH' =divH-e^5 ii (4) 

i=l 

(where N is the dimension of the modified system phase space), respectively. So, the dissi- 
pation in the modified system is greater than in the response one and it increases with the 
growth of the coupling strength e. 

The generalized synchronization regime arising in Q may be considered as a result of 
two cooperative processes taking place simultaneously. The first of them is the growth of the 
dissipation in the system an( l the second one is an increase of the amplitude of the external 
signal. Evidently, both processes are correlated with each other by means of parameter s and 
can not be realized in the coupled oscillator system independently. Nevertheless, it is clear, 
that an increase of the dissipation in the modified system J5J) results in the simplification of 
its behavior and the transition from the chaotic oscillations to the periodic ones. Moreover, 
if the additional dissipation is large enough the stationary fixed state may be realized in the 
modified system. On the contrary, the external chaotic force eAx(i) tends to complicate the 
behavior of the modified system and impose its own dynamics on it. Obviously, the generalized 
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Fig. 1 - The boundary of the GS onset on the (u>d, e) parameter plane for two unidirectionally coupled 
Rossler systems ||KJ. The oj-parameter of the response system has been fixed as uj r = 0.95. 



synchronization regime may not appear unless own chaotic dynamics of the modified system 
is suppressed. 

Note also, that the onset of GS is determined by the stability of the periodic regimes 
of the modified system (j2J) which does not seem to depend on the mismatch of the control 
parameters gd,r of the unidirectionally coupled oscillators. The stability of the periodical 
regimes is caused by the property of the modified system only. Therefore, the value of ecs 
is not supposed to depend greatly on the parameter mistuning( 1 ). Nevertheless, the sensitive 
dependence of the coupling strength value sgs corresponding to the onset of GS on the control 
parameter mismatch may take place as it has been discussed above (see also [21]). 

To explain the mechanisms causing such dependence of the coupling strength on the control 
parameter mismatch let us consider two unidirectionally coupled Rossler oscillators 

id = -UdVd - Zd, &r = ~u r y r - z r + e(x d - x r ), 

y d = Lo d x d + ay d , y r = u> r x r + ay r , (5) 

z d = V + z d (x d - c), z r =p + z r (x r -c), 

where e is a coupling parameter. The control parameter values have been selected by analogy 
with [21] as a — 0.15, p — 0.2, c = 10.0. The u r - parameter determining the main frequency 
of the response system has been selected as ui r = 0.95, and the analogous parameter oj d of 
the drive system has been varied in the range from 0.8 to 1.1 providing the slight mismatch 
of the interacting oscillators. 

The boundary of the GS regime in system (JjJ on the (u> d , e)-plane is shown in Fig. ^ 
The GS onset has been detected by the conditional Lyapunov exponent computation [10] and 
verified by the auxiliary system method [11]. The threshold of the GS arising appears to be 
essentially higher for the small control parameter detuning of the considered systems rather 
than for the large ones. At the same time, if the parameter mismatch (uj d — oj r ) is large enough, 



( 1 )This conclusion agrees well with numerical results of [23]. 



4 



EUROPHYSICS LETTERS 




Fig. 2 - The bifurcation diagram of the modified Rossler system ||SJ versus the parameter e* . The 
value of the parameter e* when the modified system starts demonstrating the periodic dynamics is 
shown by an arrow. 



the coupling strength value ecs corresponding to the GS regime onset is almost independent 
on the value u>d of the drive system (see Fig.^l. 

Such system behavior may be explained by means of the modified system method discussed 
in detail above. The response system (JJjJ can be reduced to the modified system 

ij m = uj r x m + ay m , (6) 

To separate the discussed above different mechanisms resulting in the GS arising the e* pa- 
rameter has been used instead of e. 

In Fig 2 the bifurcation diagram for the modified system © is shown. It is clear that 
with the increase of the e* parameter this system undergoes the transition from the chaotic 
to periodic oscillations through the inverse cascade of the period doubling. One can easily 
see, that starting from the value e* rj 0.06 the periodic oscillations take place in the modified 
system (JSJ). So, if the e* parameter value is large enough the autonomous modified system 
displays the periodic oscillations. 

To study the GS arising the non-autonomous dynamics of the modified system (0 

x ra — uj r y ra z m £ x m -\- £F(t*) : 

y m = u r x m + ay m , (7) 

should be considered, the time series Xd{t) of the drive system being used as the external force 
F(t). 

For the control parameters pointed above the Fourier spectrum of the drive system (JSJ) 
is characterized by the single well-defined frequency component f m . When varying the LJd- 
parameter the main frequency f m also changes. Therefore, the behavior of the modified 
system {7J) under the external harmonic signal F(t) = A cos(2ir f m t) may be considered as the 
first approximation of the dynamics of two unidirectionally coupled chaotic oscillators J5j. 
The amplitude A of the external harmonic signal F(t) should correspond to the drive Rossler 
system behavior. 

Obviously, when the main frequency of the response system and the frequency f m of the 
external signal are close enough with each other the frequency entrainment may take place 
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Fig. 3 - The synchronization tongue of the non-autonomous modified system under the external 
harmonic signal on the (cjd,s) plane. The area of the period-2 regime (marked by the letter "B") as 
well as the region of the chaotic oscillations (marked by "C") inside the synchronization tongue are 
shown by the different colors. The boundary of the GS onset (compare with Fig. has also been 
shown. 



resulting in the synchronization. In this case the behavior of the modified system (and the 
response one, respectively) may be different qualitatively inside the synchronization tongue 
and outside it. The location of the boundaries of the GS regime (see Fig. ^| confirms this 
assumption. When the parameter mismatch of the considered response and drive oscillators 
is large enough (the main frequencies are not entrained and the system @ is outside of the 
synchronization tongue) the GS regime onset is observed practically for the same value of the 
coupling strength e independently of the tod parameter value. In this case the threshold of the 
GS regime arising is defined mainly by the properties of the modified system (JJJ. Obviously, 
the features of the GS onset location for the small mistuning of the control parameters may 
be defined by the system dynamics inside the synchronization tongue. 

Let us consider the behavior of the non-autonomous modified system J7J) under the external 
harmonic signal in detail. The value of the e* parameter has been fixed as e* = 0.11 that 
corresponds approximately to the threshold of the GS arising when the control parameters 
of the coupled Rossler oscillators are detuned sufficiently. The amplitude of the external 
harmonic signal has been used in the form A = eB, the value of B parameter has been 
selected as B » 0.105 for the levels of the energy of the harmonic signal F(t) — A cos(2n f m t) 
and the main spectral component of the Fourier spectrum of the drive system Xd(t) to be 
equal. In this case the e parameter determines the intensity of the external influence on the 
modified system and may be considered as a quantity corresponding to the coupling strength 
of the coupled oscillators JHJ • Such normalization allows us to mark the synchronization region 
of the non-autonomous modified system on the (cud, e) directly. 

The synchronization tongue of the non-autonomous modified system Q on the {ojd',s) 
plane is shown in Fig. 13 The region A presents the set of the parameter values corresponding 
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to the periodic oscillation regime with the frequency f m of the external harmonic signal. With 
the increase of the e parameter value inside this region the system undergoes the transition 
to chaotic behavior through the cascade of the period doubling. The region C corresponds to 
the chaotic oscillations, while the area B presents the parameter values of the period doubling 
cascade. The boundary of the GS regime is also shown in Fig.|21for the cause of the dependence 
of the GS arising threshold ecs on the control parameter mismatch to be explained clearly. 
One can see easily, that the behavior of the response system differs essentially inside and 
outside the synchronization tongue of the modified system. 

Let us consider briefly the cause of such location of the GS onset. The external signal 
(chaotic time series of the drive Rossler system or the simulating harmonic signal) excites the 
proper chaotic dynamics of the modified system. Therefore, for the considered coupled Rossler 
systems the GS regime starts being destroyed with the growth (inside the synchronization 
tongue) of the e parameter. To suppress the chaotic oscillations in the modified system the 
value of the dissipation e* should be increased. As in the coupled systems J5J the coupling 
strength e plays simultaneously the role both of the dissipation and of the intensity of the 
drive system signal, the suppression of the proper chaotic dynamics of the response system 
takes place for the more larger values of the parameter e. Therefore, the threshold of the GS 
regime arising is shifted towards the large values of the coupling strength e. 

Comparing the boundary of the GS regime in the coupled Rossler systems (jSJ with the 
location of the synchronization tongue of the non-autonomous modified system Q under the 
external harmonic signal on the (u)d,s) plane it is necessary to take into account that the 
synchronization area has been obtained for the fixed value of the parameter e*. Moreover, 
the chaotic dynamics of the drive system influencing on the location of the GS onset has also 
been eliminated from the consideration. Therefore, this synchronization tongue on the (u;^, e) 
plane has only the qualitative character allowing to explain the mechanisms resulting in the 
GS onset features. 

In conclusion, we have explained the peculiarity of the GS onset in the unidirectionally 
coupled Rossler oscillators by means of the modified system dynamics consideration. The 
character of the GS onset location is determined by the entrainment (or, on the contrary, by 
the absence of it) of the main frequencies of the Fourier spectra of the interacting systems. 
Inside the synchronization tongue the GS onset may be shifted towards the large values of the 
coupling strength if the external influence excites the proper chaotic dynamics inside the area 
of the main frequencies entrainment. It is important to note that this conclusion has been 
made for the chaotic systems with a distinct main frequency in the power spectrum, so, the 
discussed phenomena and the presented arguments are specific to such systems. Indeed, the 
mechanism discussed above is determined by the interaction of the main spectral components 
of coupled chaotic oscillators. Probably, one can expect that the phenomenon of the increase 
of the coupling strength ecs corresponding to the onset of the generalized synchronization 
regime for the identical parameter values in comparison with the detuned ones may not be 
observed for the system with broad power spectra without distinct main frequency. 
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